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Abstract 

We present an inverse scattering approach to defects in classical integrable field 
theories. Integrability is proved systematically by constructing the generating func- 
tion of the infinite set of modified integrals of motion. The contribution of the 
defect to all orders is explicitely identified in terms of a defect matrix. The underly- 
ing geometric picture is that those defects correspond to Backhand transformations 
localized at a given point. A classification of defect matrices as well as the corre- 
sponding defect conditions is performed. The method is applied to a collection of 
well-known integrable models and previous results are recovered (and extended) di- 
rectly as special cases. Finally a brief discussion of the classical r-matrix approach 
in this context shows the relation to inhomogeneous lattice models and the need to 
resort to lattice regularizations of integrable field theories with defects. 
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Introduction 



The topic of defects, or impurities, in integrable systems has quite a rich literature, es- 
pecially for quantum aspects [H [2j El SI El EJ, [HI E] , even if quite a lot remains to be 
done. Strangely enough, the problem of integrable defects in classical field theories had 
received less attention. The pioneering paper [10] is worth mentioning for the introduc- 
tion of a so-called "spin impurity" in the nonlinear Schrodinger equation as a first step to 
tackle the problem on the half-line with integrable boundary conditions. This topic has 
been revived recently in the series of papers by P. Bowcock, E. Corrigan and C. Zambon 
[TT| [T2] . The lagrangian formalism is used in all these papers to described integrable field 
theories with internal boundary conditions interpreted as the presence of a defect. The 
defect conditions emerge from a local lagrangian density concentrated at some fixed point 
and are obtained from a variational argument. The question that is addressed is then: 
how to select conditions which leave the full theory integrable? The common underlying 
philosophy is to impose that a modified momentum, taking into account the presence of 
the defect, should be a conserved quantity while the breaking of translation invariance 
obviously entails that the bulk momentum will not be conserved. It turns out that this 
idea allows to pick up certain classes of defect lagrangian. Then, a general argument 
for integrability is based on the construction of modified Lax pair involving a limiting 
procedure. It is checked explicitely for a few conserved charges of certain models. One 
must note the nice observation made in each case: this procedure yields frozen Backlund 
transformations as the defect conditions for the fields. 

The object of this paper is to unify the results obtained by this case by case approach. 
We take advantage of the common features that have been observed. To do so, we use 
the efficient inverse scattering method formalism (instead of the lagrangian formalism) 
and implement defect conditions corresponding to frozen Backlund transformations. It 
is important to note that the role of Backlund transformations as a means to generate 
integrable boundary-initial value systems solvable by inverse scattering method has been 
discovered and used in [131 HH US]- The idea was to fold two copies of the original 
integrable system related by Backlund transformations by using compatible reductions 
on the fields (for example u{x) = u(—x)). Here, we do not fold and use the fact that 
Backlund transformations have a very nice formulation in the inverse scattering method. 
They can be encoded in matrices, representing gauge transformations of the underlying 
auxiliary problem and, in the present context, giving rise to defect matrices. Thanks 
to this formulation, we are able to prove systematically the existence of an infinite set 
of modified conservation laws, ensuring integrability. The main result of this paper is 
the explicit identification of the generating function of the defect contributions at all 
orders, i.e. for any conserved charge, for any integrable evolution equation of the AKNS 
[T6] or KN [IT] schemes of the inverse scattering method. This provides an efficient 
algorithm to compute the modified conserved quantities, given the defect matrices. One 
of the advantage of the method is that the proof of integrability does not require any 
modification of the usual Lax pair formulation for integrable field theories. Another is 
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that there is no guess work for finding the defect contributions. They are obtained from 
a classification of defect matrices. 

The paper is organized as follows. In Section 1, the general auxiliary problem formal- 
ism we use is presented. We establish our main results about the infinite set of conservation 
laws in the presence of a defect. The generalization to several defects is also explained. 
In Section 2, the defect matrices are classified within a certain class of gauge transfor- 
mations of the auxiliary problem. In section 3, we illustrate our systematic method on 
several well-known examples of integrable nonlinear equations. They correspond to all 
the classical field theories that have been explored in the lagrangian formalism (with the 
exception of the affine Toda field theories). For these models, all the previous results 
are recovered (and even generalized) and are extended to higher orders. Section 4 is de- 
voted to the extension of the method to another inverse scattering method scheme, the 
Kaup-Newell scheme [17] , which describes other classes of integrable nonlinear equations 
including e.g. the derivative nonlinear Schrodinger equation. In section 5, we discuss in 
more detail the question of integrability of such models with defects. It is argued that our 
approach allows to make a connection between the lagrangian approach and the classical 
r-matrix formalism. This requires the use of lattice regularizations. Our conclusions and 
perspectives for future investigations are gathered in the last section. 

1 General settings and results 
1.1 Lax pair formulation 

In the AKNS scheme [16] , an integrable evolution equation on the line can be formulated as 
a compatibility condition, or zero curvature condition, of a linear differential problem for 
an auxiliary wavefunction \l/(x, t, A) involving two 2x2 matrix- valued functions U(x, t, A) 
and V (x, t, A) such that 



where the subcripts x and t denote differentiation with respect to these variables. In 
the rest of the paper, we will drop the arguments whenever this is not misleading. The 
parameter A is called spectral parameter. Then, for appropriate choices of U and V, the 
integrable system at hand is equivalent to the compatibility condition ^f xt = ty tx giving 
rise to the so-called zero curvature condition 



Large classes of integrable nonlinear evolution equations can be described this way among 
which some of the most famous are the cubic nonlinear Schrodinger (NLS), sine/sinh- 
Gordon (sG), Liouville, Korteweg-de Vries (KdV) or its modified version (mKdV). It 





VA , U t -V x + [U,V}=0. 
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is known that this presentation allows one to construct generically the infinite set of 
conservation laws associated to the integrable equation. For self-containedness, we recall 
the main ideas. Let us fix U and V to be 2 x 2 traceless matrices as follows 

^(~; A i)-^'- v =( A c -a)- ™ 

where q(x,t) and r(x,t) are the fields satisfying the evolution equation. In this paper, 
we will fix the class of solutions to be that of sufficiently smooth [18] decaying fieldqj as 
| a; | — > oo and the following behaviour is assumecH 

A(x,t, A) -»• u(X) , B(x,t, A), C(x,t, A) -> as \x\ -> oo . (1.4) 

The vector-valued function \I/ is split as 

-CO- (i - 5) 

Let T = $2^1 1 j then the identification of the infinite set of conservation laws follows from 
a conservation equation 

(qT) t = (BT + A) X , (1.6) 

and a Ricatti equation for T 

T x = 2i\T + r-qT 2 . (1.7) 

This equation follows directly from the x part of ( 11.11) . The conservation equation is 
obtained from the t part to get T t , from the (12) element of (11 .2D to get qt and combining 
with (11. 7p and with the (11) element of ( 11.21) giving A x . 
Thus, expanding r as A — > oo 

n=l v ' 

the conserved quantities read 

/■oo 

I n = qY n dx , n > 1 , (1.9) 



where 

n-l 

Ti = -r , r n+ i = v nx + ^ ^ r fc r„_ fc , n>l. (1.10) 



fc=l 



-"^The decaying properties are chosen so as to ensure certain analytic properties of the scattering data 
for (jl.ip . see e.g. [TB]- Typically a polynomial decay is sufficient. 

2 The only exception is the Liouville equation for which no specific boundary condition is assumed. 
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1.2 Implementing defect boundary conditions 

Generally speaking, a defect in (1 + l)-dimensional integrable field theories can be viewed 
as internal boundary conditions on the field and its time and space derivatives at a given 
point on the line. In other words, one wants to glu together two solutions of the evolution 
equation in a specific way and at a particular point. To this end, let us consider another 
copy of the auxiliary problem. We introduce another Lax pair U, V defined as in (11.31) 
and (11.4j) with q, r replaced by q, r. We consider the analog of (11. ip for 



V(x,t,\) = L(x,t,\)V(x,t,\) . (1.11) 

The matrix valued function L(x,t,X) satisfies the following partial differential equations 
for any x and t, 

L x = UL-LU, (1.12) 
L t = VL-LV. (1.13) 

In this paper, we want to think of this matrix as generating the defect conditions at a 
specific point, xq say. Following the terminology of [HI [12], L is called the defect matrix. 
We present a classification of the simplest nontrivial such matrices in the next section. 

Now we turn to the general contruction of the generating function of the infinite set 
of modified conservation laws due to the presence of a defect. This is the main result of 
this paper and establishes integrability for any nonlinear integrable equation of the AKNS 
scheme with a defect realizing a frozen Biicklund transformation. For particular models 
(NLS, sG, Liouville, KdV and mKdV), it proves to all orders the results of [III [12] about 
the defect contribution and gives an explicit form for it. In addition, this is done without 
resorting to a modified Lax pair formalism involving a complicated limiting procedure 
to construct the conserved charges. To illustrate this, we will discuss those particular 
examples in Section [31 

To fix ideas, we choose a point xo 6 R and we suppose that the auxiliary problem 
(II. ip exists for x > xq while the one for U and V exists for x < xq. We also assume that 
the two systems are connected by the relations (jl,12p and (I1.13P at x = xq. Then, the 
following holds 

Proposition 1.1 The generating function for the integral of motions reads 

'(A) = 4:S(A) + Cf(A) + W(A), (1.14) 

where 

42(A) = r°qfdx, (1.15) 



■oo 

00 



C?(A) = / QTdx, (1.16) 
Idefect(X) = -hi(Zai + L 12 T)\ X=X0 , (1.17) 
and Lij 's are the entries of the defect matrix L. 
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Proof: From the general result (11. 6p . we get 



(qT) t = (BT + A) x , Wx>x , (1.18) 
qf) = (bT + A) , \/x<x , (1.19) 



where T is defined from \l/ as in the previous section, T = ^ 2 ^i ■ From this and the 
rapid decay of the fields, we get 

f°° f x ^ / _ v 

8 t / qTdx + d t qTdx = [BT + A- (BY + A)) \ X=XQ (1.20) 

J XQ J— OO 

The crucial point now is that the right-hand-side is a total time derivative of a quantity 
evaluated at x = xo: it is the contribution of the defect to the conserved quantities as we 
now show. From (11.111) we get T = (L 2 i + -^22r)(Ln + Li 2 r) _1 . Then, using (11. 13ft at 
x = xq to eliminate A and B, one gets 



BT + A- (BT + A)) \ X=X0 (1.21) 
= {d t L n + d t L 12 T + L 12 (C - 2AT - BY 2 )} (L n + L 12 Vy l \ x=X0 . (1.22) 

The final step consists in noting that the t part of (11.11) implies another Ricatti equation 

T t = C-2AT-BT 2 , (1.23) 

so that 



Xq J —OO 



qTdx = ^ 11 + Ll2T h x=X0 , (1.24) 



Therefore, 

d t I(X) = 0. (1.25) 



Note that this holds in all generality for any integrable evolution equation in the AKNS 
scheme with decaying fields. The latter can be relaxed for some models, and we discuss 
in detail below the Liouville equation for which this does not hold. 



1.3 Several defects 

It is quite straightforward to repeat the general argument for the construction of conser- 
vations laws in the case of several defects. Suppose we have iV + 1 auxiliary problems 
connected two by two at points x±, . . . , x^- by matrices L 1 , . . . , L N . Then, it is easy to see 
that the contribution of these N defects is the sum of the contributions from each defect. 
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Indeed, using xq = — oo, xn+i = +00 and otherwise obvious notations, the generating 
function for the integral of motion reads 

N+l N 

/(A) = E 7 W A )+£WA), (1.26) 

i=i i=i 

HmW = r U^dx , j = l,...,JV+l, (1.27) 

4/«*(*) = ln((L i ) 11 + (L J ) 12 r i )U =:co , j = l,...,iV. (1.28) 



2 Defect matrices 

In this section, we derive a large class of defect matrices satisfying (j!.12p and (11.131) 
together with the associated conditions they entail on the fields: the Backhand transfor- 
mations. The latter will become the defect conditions when imposed at x = xq. 



2.1 Generalities 

The matrix L preserves the zero curvature condition as is easily seen by writing L xt = L tx , 

(VA, U t — V x + [U, V] = 0) & (VA, U t - V x + [U, V] = 0) . (2.1) 

In other words, if q, r are solutions of the evolution equation described by (U, V) then 
q, r are solutions of the evolution equation described by (U, V) under the transformation 
induced by L and vice versa. This is just the usual definition of a Backhand transformation 
and this shows the connection with the idea of frozen Backhand transformations discussed 
above. In other words, we look for defect matrices in the class of matrices realizing 
Backlund transformations between two nonlinear integrable evolution equations. Note 
that the evolution equations need not be the same in general. If they are, the terminology 
auto-Backlund transformation is usually used. Such matrices are sometimes referred to 
as Darboux matrices (see e.g. [IH]). Even if a lot is known on these matrices, we proceed 
with their derivation in the form needed for this paper. We adopt a pedestrian method 
which does not require any previous knowledge of their theory. In particular, no reference 
to the wavefunction of the auxiliary problem or to a special Riemann problem is needed 
(which are usually the methods encountered in the literature). 

Let us establish some general facts about L. First, there is some freedom in its normal- 
ization coming from the invariance of the zero-curvature condition under the transforma- 
tion (U,V) —> (M~ l UM, M~ l VM) for any invertible matrix M independent of x and t. 
This also obviously preserves the tracelessness property. In particular, left multiplication 
of L by M~ l amounts to apply this transformation to (U, V) while right multiplication 
by M applies it to (U, V). Then, we have the 
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Proposition 2.1 The determinant of L is independent of x and t 

det L(x,t,X) — f(X) . (2.2) 

Proof: The result follows from the Jacobi formula 

(detL(x,t,X)) x = det L(x, t, X)Tr(U - U) , (2.3) 
(detL(x,t,X)) t = det L(x, t, X)Tr(V - V) , (2.4) 

and the tracelessness of U, U, V, V. ■ 

At this stage, it is hard to go further without specifying U, U, V, V a bit more. Let us 
simply note that given U, U, V, V and initial-boundary values for the fields, the integration 
of (11 .121) and (11.131) gives for instance (the path of integration being irrelevant) 

L(x,t,X) = L(x ,t ,X)+ [ X (UL-LU)\ T=t dy+ [ (V L - LV)\ y=Xo dr . (2.5) 

Jx J t 

The formal iteration of the previous equation suggests that, in general, L has a compli- 
cated Laurent series structure as a function of A. In the following, we will assume that L 
has only a finite number of terms and, recalling that it is defined up to a scalar function 
in A, we will look for a solution of the form 

N 

L{x,t,X) = J2 L(n) ( x ^) X ~ n ■ ( 2 - 6 ) 

n=0 

Actually, we shall consider the case N = 1 which we study in detail. For convenience, we 
also restrict our attention to auto-Backlund matrices. We comment later on on the fact 
that this is not necessary in our approach, one of the crucial ingredient being simply that 
the evolution equations have the same dispersion relation (see (12 . 181) below). 

2.2 Construction for N=l 

The defect matrix is of the fornj§ 

L(x,t,X) = L(°\x,t) + LW(x,t)X- 1 . (2.7) 
The defining relation (11.121) is equivalent to 

= [L<Va] , (2-8) 
= ip l \a z ]+WL^ -L^W , (2.9) 
L« = WL^-L^W. (2.10) 



3 Note that it is implicitly assumed that both L( ) and are not trivial since otherwise, the Backlund 
transformation is essentially the trivial one W = W. 
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If V,V are polynomials in A with coefficients V^\x, t), V^(x, t), j = 0, • • • , N, equation 
(11.131) is equivalent to 

L (D = y(0) L (l) _ L (l)y(0) ? (2 U) 

L f (0) = y«L« - L^V^ + V<®L<® - L®V® , (2.12) 
= V^L^-L^V^ + V^L^-L^V^, (2.13) 

i (2.14) 
= VWL®-L<®VW. (2.15) 

If V,V are polynomials in A" 1 , the equations are the same under the exchange L^ ' <-> L^. 

Let us make a few remarks. First, when we have found the matrix L, equations (12.101) 
and (12.111) will give the x and t parts of the corresponding Backlund transformations for 
the fields. Then, in traditional approaches, equation (12. 9p is used to construct new soliton 
solutions, W, from given solutions W and the knowledge of the Backlund transformation. 
We will not discuss this last step in this paper and refer the reader to the vast literature 
on the subject (see e.g. [20] and references therein). 

We now proceed with the statement of the general results of this section. 

Proposition 2.2 Defect matrix 

The defect matrix L has the following general form 



L=X 2 + \-i\ ^ J i n , (2.16) 

a 3 i | a + =F [at - 4a 2 a 3 \ 



where 



a 2 = --(q-q) , a 3 = -(f - r) , (2.17) 



and a± G C are the (x,t-independent) parameters of the defect. 

Proof: Equation (12.81) implies that in (12.71) is diagonal and then, equation (12.91) shows 
that the diagonal elements do not depend on x. Therefore, we can consider equation 
(12.121) as \x\ — > oo. Recall that we have 

V,V^Lo(X)a 3 , |z|->oo. (2.18) 

N 

Writing tu(X) = V]o; (n) A n and denoting = lim L (1) (x,i), we get 

n |a:|->oo 
n=0 



L? = ^[a 3l L£]+JV[a 3 ,L^] (2.19) 

= -iu {l) lim {WL^ -L {0) W) (2.20) 

\x\— too 

= 0, (2.21) 



S 



where we have used equation (12.81) and (12.91) in the second equality and the fact that we 
consider decaying fields in the last equality. The proof is similar if V, V are polynomials 
in A -1 . So, as explained above, we can left multiply with (-f^ -*) and work with L' = 

(L( )) 1 L. We drop the ' in the following but remember that should now be II2 in 
all the equations (l2.8H2.15j) . Next, denote 



a 3 



a -2 
<i 1 



and cti, a 2 its eigenvalues. Then, equation (12.91) gives immediately a 2 

a 3 



r). Now, the elements a± and (14 are easily computed from 
aia4 — 0203 = a\a 2 and a\ + 04 = a\ + a 2 , 



(2.22) 

■3(9 - q), 

(2.23) 



and introducing ot± = a\±a 2 . Finally, we need to show that a± is independent of x and t. 
Let i\ and £ 2 be the eigenvalues of L. It is enough to prove that i\ and l 2 are independent 
of x and t. In turn, it is sufficient to prove that £\£ 2 and l\ + £ 2 are independent of x and 
t. From proposition 12. 11 we already know that £\£ 2 = /(A). Next, we prove £i + £ 2 = #(A). 
Suppose £1 + £ 2 = g(x,t, A), then, using (11.121) and (12.91) 



9x 



Tr 



L(W - W) 



A _1 Tr 



L {1) (W - W) 
{Tr [L« [(73, LW]] 



z 

A 
0. 



Now, (2 t = Tr 
So & = 0. 



L(V - V) 



can be evaluated as x — > 00 for which we know that V- 



(2 


24) 


(2 


25) 


(2 


26) 


(2 


27) 


V - 


■* 0. 




■ 



The Backhand transformations associated with the matrix L read: 

• For the x part, 

aix = 9«3 - ra 2 , 
a 2x = qa 4 - qa\ , 
a 3x = rai - ra 4 , 
a 4x = ra 2 - qa 3 , 

• For the t part if V, V are polynomials in A, 

a« = (A^-A^)a 1 + B^a 3 -C^a 2 , 

a 2t = (A^+A^)a 2 + B^a 4 -B^ ai , 

a 3t = -(^+A^)a 3 + C^a l -C^a A , 

a u = -(A^-A^)a 4 + C^a 2 -B^a 3 , 



(2.28) 
(2.29) 
(2.30) 
(2.31) 

(2.32) 
(2.33) 
(2.34) 
(2.35) 
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For the t part if V, V are polynomials in A 1 , 

2A (0) a 2 + 5 (0) (a 4 - a x ) + d 2 5 (1) - di£ (1) , 
a 3 + C (0) (ai - a 4 ) + diC (1) - d 2 C (1) 



ait 

03t 



«4t 



«3 



(2.36) 
(2.37) 
(2.38) 
(2.39) 



where 01,02,03,04 are as in Proposition 12.21 At this stage, it seems that there is an 
overdetermination since there are four equations for each part whereas only two of each 
type are needed (the x and t transforms relating q and q and those relating f and r). It 
turns out that half of them are indeed redundant. 

Proposition 2.3 Backlund transformations 

The Backlund transformations corresponding to L are given by the equations for a 2 and 
a 3 in $2.2<k\2.30\) and \2.33\\2.3$ or \2.3T\\2.38\) . The remaining equations for a\ and a 4 
can be deduced from them. 



Proof: We know that the eigenvalues of are constant so 

(oia 4 )x = (a2fl3)x- , a lx + a Ax = . 
From this and the equations f l2.29ll2~30l) for a 2 and 03, we deduce 
ai x (a A - ai) = a A (qa 3 - ra 2 ) + a 1 (fa 2 - qa 3 ) • 
Now, using (q — q)a 3 + (f — r)a 2 = 0, 

(a ix - qa 3 + ra 2 )(a 4 - a x ) = . 



(2.40) 



(2.41) 



(2.42) 



The possibility a 4 = a 4 must be rejected in general since together with a± x + a^ x = it 
would imply that ai and a 4 are independent of x. Thus, we obtain the equation for ai 
and hence for a 4 . 

The proof for the t part is similar. Useful identities in getting the result are obtained 
from 



Tr L t = Tr 
and expanding in powers of A or A -1 . 



L(V - V) 



0, 



(2.43) 



We finish this general discussion by making a connection with Darboux matrices (see 
e.g. [IE]). Suppose that a± 7^ a 2 then we can define 



a 2 — ai 



■(L«-ailI 2 ) 



(2.44) 
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and multiply L(X) by (since it is denned up to a function of A) to get 

L(A) = JI 2 + ^^P. (2.45) 

The important fact is that P is a projector (P 2 = P can be checked directly from the defi- 
nition (I2.44p ). The form (12.451) for L is usually encountered where the so-called Backlund- 
Darboux transformations are used to generate multi-soliton solutions from a given one. 
This form is also useful to exhibit the inverse of the Backlund matrix 

L -i(A) = i a -£lZ^p. (2.46) 

A + Oi2 



3 Examples 

In this section, our results are applied on a variety of examples. We are able to reproduce 
the results of [TTl [T2] in a very simple way. For the NLS equation, we obtain a more 
general result corresponding to a Backlund transformation with two real parameters (as 
it should, see e.g. [21]) instead of one. For the Korteweg-de Vries (KdV) and modified 
KdV equations, we obtain the defect contribution directly in terms of the original fields 
and reproduce the lagrangian approach expressions in terms of "potential" fields. For each 
model, the defect conditions are given and consistently reproduces the associated well- 
known Backlund transformations, but taken at x = xq here, as expected by construction. 
We gather the examples in three classes according to certain symmetry considerations 
yielding information on the defect parameters a±. 



3.1 Class I: q = u, r = eu*, e = ±1, u complex scalar field 

Let us introduce 



K 



1 

e 



We have the following symmetries (we drop x and t) 

U*(X*) = KU^K- 1 , *7*(A*) = KUiXjR- 1 , 



(3.1) 



(3.2) 



and we assume that V and V have the same properties. Therefore, we can look for the 
Backlund matrix L such that L*(A*) = KL(X)K~ l . This implies «2 = ct* so a + 6 K. and 
a_ G iR. Then, it can be shown that the remaining conditions imply that the nontrivial 
Backlund matrix reads 



L(A) = fl 2 + A" 1 



( \ |a+ ± iy/P 2 + e\u - u\ 2 \ 



- u) 



era — it 



, (3.3) 
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where (3 = ict- G M. For this class, the defect matrix, and hence the defect conditions 
are parametrized by two arbitrary real numbers a + and j3. In the case e = — 1, we see 
that the transformation is such that \u — u\ 2 < (3 2 . We must take into account the fact 
that r = eg* in the discussion of the integrals of motion. In particular, it turns out that 
to generate real integral of the motion (real classical observables) , one has to consider the 
following combination 

r ym (X) = - J*(A*)) . (3.4) 

So in practice, we will compute the contribution of the defect as 

QJectW = ~i (HLn + L l2 T) - ln(L n + L 12 r)*) \ X=XQ . (3.5) 

Now we expand (13.41) in powers of A -1 up to order 3 to illustrate the method. For 
convenience, we define il e = \J ft 2 + e\u — u\ 2 . At order A" 1 , we find that the modified 
conserved density reads 

Xo POO 

\u\ 2 dx + / \u\ 2 dx eil e \ x=xo , (3.6) 



where the last term is the explicit defect contribution. Similarly, at order A -2 , the modified 
conserved momentum is 

/xo roo 
i(uu x — u*u x )dx + / i(uu* — u*u x )dx — i [{u*u — uu*) =f 2e/3f2 e ] \ x=xo . (3.7) 
-oo J xo 

Finally, the modified conserved energy is found to be 

xo roo 

{\^x\ 2 + e\u\ A )dx + / (\u x \ 2 + e\u\ 4 )dx 

oo J xo 

+ [^a(N 2 + \u\ 2 ) =F ^a(3/3 2 - n 2 e ) - i{3(uu* - uu*)j \ X=X0 . (3.8) 

These results hold for any member of class I so in particular they hold for the cubic 
focusing (e = —1) or defocusing (e = 1) nonlinear Schrodinger equation for the complex 
scalar field u 

iu t + u xx = e\u\ 2 u , (3.9) 
and similarly for u. Indeed, this equation is obtained in the AKNS scheme by taking 

A(X) = -2i\ 2 + i\u\ 2 , B(\)=eC*(\*) = 2\ + iu x , (3.10) 
and similarly for u. 

Now the corresponding defect conditions can be derived from the general Backlund 
transformations given in the previous section. From the symmetry 03 = ea*,, we need only 
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consider (12.291) and (I2.33[) . Note that the x part of the defect conditions is the same for 
all the models in class I. Here, for NLS, we have at x = xq 



{u — u) x = ia + (u — u) ± (it + u) ^/ j3 2 + e\u — u\ 2 , (3-H) 
(u — u) t = —a + {u — u) x ± i(u + u) x \J (3 2 + e\u — u\ 2 + i(u — u)(\u\ 2 + \u\ 2 ) (3.12) 

Setting e = — 1 and choosing the — sign in (13.61) . (13.71) . (13.81) . we note that we have to 
impose further (3 = to recover the results of [12] (where the notation f2 = a/« 2 — \u — u\ 2 
is used). This is due to the fact that the lagrangian the authors took for the defect (the 
E functional in their notations) is not the most general one. It corresponds to particular 
Backlund transformations with (3 = 0. 

3.2 Class II: q = u, r = ew, , e = ±1, u real scalar field 

This class is a subclass of class I with u* = u (and u* = u). this immediately implies 

a+ = 0. (3.13) 

The defect matrix for this class reads 



L(X) = 2 2 + A- ( ^ 2 + *(« " -) 2 ) . (3 . 14 ) 

V \e(u-u) =Ff Va 2 + e(u - u) 2 J 

So we can simply use the result of the previous class, setting u* = u and u* = u (and thus, 
without symetrizing) . We exhibit the first orders for a few examples, taking advantage of 
specific forms of the defect matrix in each case. 

3.2.1 Modified Korteweg-de Vries equation 

The modified Korteweg-de Vries equation equation 

u t - 6eu 2 u x + u xxx = , (3.15) 
is obtained in the AKNS scheme by taking 

A(X) = -4zA 3 - 2i\eu 2 , B(\) = eC*(A*) = 4X 2 u + 2i\u x — u xx + 2eu 3 , (3.16) 
and similarly for u. The modified conserved density reads 

x 



u 2 dx+ I u 2 dx =F ey a 2 + e(u — u) 2 \ x=X0 . (3-17) 



The modified conserved momentum is 

rXQ poo 



/x poo ^ 
uu x dx + / uu x dx (u 2 — u 2 + a 2 )\ x=XQ , (3.18) 
OO iln ^ 



-OO J Xq 
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as can be directly checked by integration by parts (the constant a 2 being irrelevant). 
Finally, the next order yields 

/XQ r-co el 

{u 2 x + eu 4 )dx + j (u 2 x + eu A )dxT^e (u 2 + u 2 ) - \ x=Xo , (3.19) 
-oo J xo 3 J 

where here f2 e = \/ a 2 + e(u — u) 2 . The corresponding defect conditions read 

(u-u) x = ±(u + u)y/a 2 + e(u - u) 2 , (3.20) 
{u-u) t = ± {2e(n 3 + u 3 ) - (u + u) xx ) ^a 2 + e(u - u) 2 . (3.21) 

It is worth noting that everything is expressed directly in terms of the fields u and u. This 
should be compared with the Lagrangian approach of [12] where this was not possible. 
The use of "potential" fields p and q such that u = p x and u = —q x is required in this 
formulation. Under this substitution, an alternative form of the defect matrix can be 
derived 

L(X) = n 2 ± A" 1 — ( C ° s( r V \ - Sil f„ - V \) . (3.22) 
2 V — sm(i; — v) — cos{v — v) J 

and we consistently recover their result for the defect contribution to the first conserved 
quantity (setting e = —1), that is e^J a 2 + e(u — u) 2 \ x=xo becomes 

a(cos(p -q)-l) \ x=xo , (3.23) 

(a constant can always be added). 

3.2.2 Sine/sinh-Gordon equation 

This example illustrates the case V, V polynomials in A -1 . The sine-Gordon equation in 
light-cone coordinates 

v xt = sinv, (3.24) 

is obtained by setting u = — ^ e = —1 and taking 

. . , , i cos v „ , , . i sin v . 
A(X) = — , B(X) = — , (3.25) 

and similarly for v. For this model, the defect matrix takes the nice following form 

ia ( cos^ -ffln- B+u 



L(A) = I 2 ± — . I +V _ ,% v , (3.26) 



2A V ~ sin 2f2 - cos 2 
where a is a nonzero real parameter. Therefore, the modified conserved momentum reads 
1 f x ° _ 1 f°° v + v 

- / v 2 dx-\ — / v 2 .dx±acos \ x =xn ■ (3.27) 

4 loo 4 J XQ x 2 ]x Xo K J 
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This is just the result in [30] for the momentum and energy but expressed here in light-cone 
coordinates. At the next order, as can be anticipated from (I3.18p . the bulk contribution 
can be integrated by parts and combines nicely with the defect contribution, leaving the 

2 

constant — ^- as the conserved quantity. Finally, the third order conserved quantity is 



1 r\ 2 V* , 1 f°° , 2 V* a V + V f 2 2 2m 

2 J (-Kx + f)dx + -J (-v 2 xx + ^)dx±-cos^-(v 2 x + v x v x + v 2 x + 2a 2 )\ x=x i3.2l 
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" J — oo * " •/ xo 

The defect conditions at x = x are given by 

(v-v) x = ±2asin^^, (3.29) 

+ = ±-sin^^. (3.30) 
a z 

3.2.3 Liouville equation 

This is another example of evolution equation obtained from V and V polynomials in 
A -1 . We need to discuss this equation in detail since, as mentioned above, the condition of 
vanishing fields at infinity and the boundary conditions (11.41) are not applicable. However, 
it is straightforward to prove that Proposition 12.21 is still valid. Also, it is still possible to 
prove Proposition II . II with a slight modification, as we now show. The Liouville equation 
in light-cone coordinates for the field v 

v xt = 2e v , (3.31) 

is obtained by taking 

u= V f , 6=1 , A(A) = ^ = -B(X) , (3.32) 

and similarly for v. The conservation laws (11.18P and (11.191) still hold since they do not 
depend on the boundary conditions. But now, (11.201) becomes 

d t / qTdx + d t qTdx = Urn (BY + A) - lim (BY + A) (3.33) 

/ / x— >oo X— >— oo 

J xo J — oo 

+ (Sr + A- (BY + A)) | x . =xo . (3.34) 

The last term in the right-hand-side is treated as before. The point is to recast the other 
two terms as time derivatives. For this equation, one has A = C = —B so BY + A = 
A(l — Y) and the Ricatti equation (I1.23P becomes 

Y t = A(l-Y) 2 . (3.35) 

Therefore, since r / 1 (this can be seen from (11.71) ) 

BY + A = -fifcln(l-r). (3.36) 
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The same result holds for T. Therefore, Proposition 11.11 is modified to the following. 
The generating function for the integral of motions of the Liouville equation reads 

'(A) = S(A) + G h< (A) + /, e/erf (A), (3.37) 

where 

v x Tdx - lim ln(l - f ) , (3.38) 



4«{fc(A) - 


1 

2 


rx 
J — oo 


T ri 9 ht (\\ — 

hulk \ A ) - 


1 

2 


POO 

J 


IdefectW — 




ln(L 



x— >— oo 



v x Tdx + lim ln(l - T) , (3.39) 

x— >oo 

ii + Li 2 r)| x=X0 , (3.40) 

and Ljj's are the entries of the Backlund matrix L. The additional contributions essentially 
kill terms arising from trivial integration by parts in the integrals of motion. 
The defect matrix can be written 

i'j ± v±v ( 1 — 1 



L(A) = H 2 + J ^* ^ j, (3.41) 

where 7 is a nonzero real constant. At first order, the modified conserved momentum 
reads 

1 r x ° 1 f°° v+v 

- / ^.dx + v x \-oc + x / vldx-v x \ oo -je ±E 2 E \ x=x0 . (3.42) 

^ ./-oo ^ Jxo 

As is now customary, the next order combines nicely to produce 

(v 2 x - 2v xx )\ 00 - (v 2 x - 2v xx )\- <x , (3.43) 
which can be checked directly to be a constant. Finally, at the third order, we have 

1 r° /_ 2 v x \ J / i)3 _ _ 

u** + — ax - ?w - — - v x u 



2 J-oo V 4 / V 6 

1 r 00 / 2 / vi 

+ 2 J \ xx + ~AJ + V X:r:r ~~ "6~ _ VxVxX 
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The defect conditions are given by 



7 e ±v -¥ (v 2 x + v 2 x + v x v x ) \ x=x0 . (3.44) 



v-v) x = 7e ± 2 , (3.45) 

(v + v) t = ±-e^(e ii -e v ), (3.46) 
7 
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3.3 Class III: q = u, r = e, , e = ±1, u real scalar field 

For this class, there is no special symmetry. The derivation of the defect matrix deserved 
special attention for this class. Indeed, if we assume that (U, V) has the same form as 
(U, V), as we have done so far, then an immediate consequence is a 3 = and a\ = a 4 = ati 
are constant. From this, equation (12.281) implies u = u i.e. we get the trivial defect 
conditions. A solution to this problem is to left multiply L by 03 or equivalently to 
consider (asUa^, a^Va^). This amounts to change the sign of the off-diagonal terms. 
Then, equation ( 12.301) implies a + = 0. Finally, taking into account the reality of the field, 
one gets the defect matrix for this class 

L(X) = I 2 + A" 1 ( ±WP 2 + Hu + u) Uu + u) \ 

K J V ~ ie HV/ 3 2 + 2e(n + n) J K J 

We apply our method to the Korteweg-de Vries equation 

u t - 6euu x + u xxx = , (3.48) 
which is obtained in the AKNS scheme by taking 

A(X) = -4iA 3 - 2ieXu + eu x , (3.49) 
B(X) = 4X 2 u + 2i\u x + 2eu 2 - u xx , (3.50) 
(7(A) = 4eX 2 + 8u, (3.51) 

and similarly for u, with appropriate change of signs. We make direct use of (j!.14p . The 
first nontrivial order is A -3 and for the modified conserved density reads 

-/ u 2 dx+- u 2 dx =f - y/f3 2 + 2e(u + u) (e(u + u) - (3 2 ) \ x=Xo . (3.52) 
2 J -00 2 < / a . 6 

The defect conditions are given by 

(u + u) x = ±{u-u)^(3 2 + 2e(u + u) , (3.53) 

(« + u) t = ± (3e({i 2 - u 2 ) -{u- u) xx ) V/S 2 + 2e(u + u) . (3.54) 

Once again, we note that everything is expressed directly in terms of the initial fields 
u and u while this was not possible in the lagrangian approach. It is easy to make contact 
with the latter by setting u = q x and u = p x . Then, 

(p + q), = ^[(p-q) 2 -P 2 } ■ (3-55) 

Taking e = 1 and setting (3 2 = —4a, we recover the result of [12] for the defect contribution 
to the density (the momentum in their setting) 

2 J Pl dx + 2 / q * dx ' ( 3 - 56 ) 

that is 
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3.4 Remarks 



We will not go into the analysis of the higher N case in (12. 6p . We simply note that a large 
class of higher N defect matrices is provided by products of N = 1 defect matrices. The 
corresponding defect conditions are then simply compositions of the defect conditions 
we derived above. In other words, the defect matrices we constructed have a group 
structure, as can be seen directly from (11.121) and (11.131) . Furthermore, Bianchi's theorem 
of permutability (see e.g. [19]) implies that this is an abelian group. An important 
question concerns the existence and properties of those higher N defect matrices which 
do not factorize as products of N = 1 ones. Their study would shed new light on possible 
new defect conditions for the well-known systems we discussed. Also, the question of 
defect matrices preserving integrability but which do not fall at all in the class discussed 
here remains entirely open. In this sense, no claim of uniqueness of defect matrices is 
made and one should remember that the proposed approach here is sufficient to ensure 
integrability in the presence of a defect. The issue of finding necessary defect conditions 
for integrability is not answered. It should also be noted that following the linearization 
argument of [HI [12], it appears that the defect conditions constructed here allow for pure 
transmission only. This has been checked explicitely for all the examples given above. 

4 Extension to another scheme 

The Kaup-Newell (KN) [17] scheme goes along the same steps as the AKNS scheme to 
produce integrable evolution equations with the essential difference that the matrix U 
involved in the x part of the auxiliary problem ( 11.11) has the following form 



One well-known model obtained in this scheme is the derivative nonlinear Schrodinger 
equation 



It turns out that our method works for this scheme too with the appropriate modifications. 
We proceed as before by taking two copies of the auxiliary problem related by a matrix 
M such that 




(4.1) 




(4.2) 



(4.3) 



M x 
M t 



U KN M - MU KN , 
V KN M - MV KN . 



(4.4) 
(4.5) 



Then we assume that the two copies are related by M at some point x = xq. 
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Proposition 4.1 The generating function for the integral of motions reads 

'(A) = 42(A) + iffi* (A) + I d efect(X) , (4.6) 

w/iere 

4S(A) = / AgTdx, (4.7) 

J oo 
✓•oo 

= / AgWx, (4.8) 
IdefectW = -ln(Mu + M 12 r)| (4-9) 

and Mjj 7 s are i/ie entries of the Backlund matrix M. For this scheme, T and T have a 
different expansion as A — > oo 

oo r 

r ^#r- (4-10) 

n=0 v 7 

with 

n 

r = -r , r n+1 = 2?r n:c + g r p r n _ p , (4.11) 

p=0 

and similarly for T. 

Proof: The proof is the same as that of Proposition 11.11 the only difference being that the 
conservation equations now read 

for), = \(BT + A) x , Vo;>Xo, (4.12) 
(gf) t = i(5f + 2^ , \fx<x , (4.13) 

due to the different A dependence of Ukn- The latter is also responsible for the different 
series expansions of T and T coming from the following Ricatti equation 

T x = Ar + 2t\ 2 T - XqT 2 , (4.14) 
and similarly for T. m 

To apply this to specific models, we would need some sort of classification of the 
matrices M along the lines of what is available for L. However, the author is not aware 
of such results. It is a problem for future investigation. 
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5 Discussion of integrability 



So far, we have shown how the infinite set of conservation laws is modified by the presence 
of a defect described by a matrix for any evolution equation falling into the AKNS or KN 
schemes. The role of an infinite set of conserved quantities is well-known in the construc- 
tion of action-angle variables in the inverse scattering method for evolution equations on 
the line. In turn, this construction ensure the integrability of the system in the sense that 
in terms of the new variables, the evolution in time is very easy to solve. Then, using 
the inverse part of the method (the Gelfan'd-Levitan-Marchenko equations [2"2"l |2"3"] ) one 
can deduce the time evolution of the original fields (see e.g. [IS])- So, in this sense, the 
systems with defect we have considered are integra 

A traditional complementary view is to reformulate the evolution equation as Hamil- 
tonian systems with a Poisson structure. In this formalism, the idea is to show that the 
conserved quantities previously constructed form a commutative Poisson algebra contain- 
ing the Hamiltonian which generates the time evolution of the fields. Then, one talks 
about integrability in the sense of Liouville. One of the advantages of this reformulation 
is the possibility of quantization and this was the object of the quantum inverse scattering 
method, see e.g. [21]. 

The method of the classical r-matrix [25] has proved very useful and fundamental 
in the discussion of these issues for classical integrable systems on the whole line. The 
situation on the half-line is also well understood [26] . Our aim in this section is to discuss 
the situation with a defect. 

Let us first recall some facts about the method of the classical r-matrix for ultralocal 
models. The basic ingredient is the 2x2 transition matrix T(x,y,\), x < y, defined 
as the fundamental solution of the x-part of the auxiliary problem at a given time (not 
explicitely displayed in T) 

d y T(x,y,X) = U(y,t,X)T(x,y,X) , T(x, x, A) = H 2 , (5.1) 

where U is given as in (11.31) . The important result reads (see e.g. [27] ) 

{Ti(x,y,X),T 2 (x,y,fi)} = [r 12 (A - fi),T 1 (x,y,X) T 2 (x,y,/J,)) , (5.2) 

where r 12(A), the classical r-matrix, is a 4x4 antisymmetric solution of the classical Yang- 
Baxter equation [2H]. For systems on the circle, it allows to show that the coefficients 
of the series expansion of TrT(x, y, A) in A are in involution and the Hamiltonian is one 
of the them. For systems on the line, the principle is the same but one has to take the 
infinite volume limit of (15.21) . usually with special care. 

For the problem with defect, we would like to mimic this procedure. The main dif- 
ference here is that there is something nontrivial going on at a single point x = Xq and 

4 For completeness, one should perform the inverse scattering method and identify the action-angle 
variables in this context. This is left for future work. 
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characterized by L (or The standard procedure becomes ill-defined in the continu- 

ous case as it involves singular Poisson brackets of the type "5(0)" where 5(x) is the Dirac 
distribution. 

Indeed, the analog of the transition matrix for x < x$ < y is 

T x °(x, y, A) = T(x Q , y, X)L~ l (x , t, X)f(x, x , A) , (5.3) 
and the task of computing 

{T?°{x,y,\),T^{x,y,ii)}, (5.4) 

involves computing {L]" 1 (xo, t, A), (xo, t, A)}. There does not seem to be a direct ap- 
proach starting from the explicit form of L as classified in this paper. However, in the 
context of finite-dimensional integrable systems, important results were obtained by E. 
Sklyanin in [25] concerning the canonicity of Backlund transformations in the formalism 
of the r-matrix approach. Transposing the results in the present context and assuming we 
have the same theories on both sides of the defect, i.e. T(x,y,X) and T(x',y',X) satisfy 
(15. 2j , for x,y > x and x',y' < x respectively, it seems reasonable to postulate that L^ 1 
is just another representation of the Poisson algebra (15.21) so that 

{L^(x , t, A), L^ixo, t,X)}= [7*12 (A - fi), L^(x , t, A) L^(x , t, A)] . (5.5) 

From this it immediately follows that 

{ y, A), T*°(x, y,n)}= [r 12 (A, fi),7?°(x, y, A) T 2 ^(x, y, fx)} . (5.6) 

Note that the generalization to defects can be described in this formalism as well. 
Using the same notations as Section 11.31 one construct the monodromy matrix 

T N (X) = T N+1 (x N ,x N+1 ,X)(L N )-\x N ,t,X)...(L 2 )-\x 1 ,t,X)T 1 (x u x 2 ,X) 

x{L 1 )- 1 (x 1 ,t : X)T°(x ,x 1: X), (5.7) 

with all the L J 's satisfying (15.51) . 

This discussion brings us to the connection between the lagrangian approach of [HI [12] 
to integrable defects and a quite standard procedure to implement inhomogeneities or 
impurities in discrete integrable systems. Indeed, through the approach of this paper, 
one can reformulate the lagrangian approach in terms of a transition matrix made of 
two bulk parts and a localised defect part realising a different representation of the same 
Poisson algebra. But this is exactly what is usually done to implement so-called impurities 
or inhomogeneities in discrete integrable systems. Let us consider a discrete ultralocal 
system of length i with N sites. The transition matrix T(A) is then a product of local 
matrices tj(X), j — 1, . . . , N satisfying 

{tjo(X),t ko ,} = S jk [r 00 '(A - n),t j0 (X)t k0 ,({i)] , (5.8) 
5 For convenience, in the rest of the paper, we stick to a single notation L for the defect matrix. 
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in the same representation (the index j represents the space of dynamical variables while 
0, 0' are auxiliary spaces, C 2 here). One can check then that T(A) = tjv(A) . . . £i(A) 
satisfies 



{r (A),T(y(/i)} = [t w (\-h),T {X)M»)] , 



(5.9) 



i.e. exactly (15.21) . To introduce an inhomogeneity at site jo say, one then chooses a 
different representation tj (X) of the same algebra. This does not change the properties of 
T(A) = t N (X) . . . tj (\) . . . ti(X) but influences the physical quantities (e.g. the integrals of 
motion) that can be computed since the latter depend on the representation at each site. 
Again, the introduction of several inhomogeneities can be done straightforwardly. This 
actually provides a solution to the above problem of computing Poisson brackets of defect 
matrices by considering lattice regularizations of integrable field theories and changing 
representations appropriately at local sites to generate defects. 

From this point of view, one can anticipate the outcome of the quantization of this 
approach. It is known that (I5.9P is the classical (H — * 0) limit of the quantum Yang-Baxter 
algebra 



and Tj(A) — > tj(X) in the h — > limit. So the quantum defect matrix C(X) encoding the 
defect conditions will satisfy the quantum Yang-Baxter algebra. This gives some support 
to the ad hoc quantization procedure adopted in [5U| for the sine-Gordon with integrable 
defect. We would like to stress that the above programme of discretization has been 
completed for the defect sine-Gordon model in the important paper [31], both at classical 
and quantum level. The approach is based on the notion of ancestor algebra |32j . 

Conclusions and outlook 

In this paper, we have reformulated the lagrangian approach to the question of integrable 
defects in the language of the inverse scattering method, taking advantage of the common 
features that had been observed on a case by case study: frozen Backlund transformations 
as defect conditions ensure integrability. The reformulation allows a systematic proof of 
this as well as an efficient computation of the modified conserved quantities to all orders 
in terms of the defect matrix. The latter, and the associated defect conditions, can be 
classified and we performed these computations for a certain class of matrices. Taking 
particular examples, we recovered and even generalized all the previous results obtained by 
the lagrangian method. It should be emphasized that this procedure provides a sufficient 



R 00 '(X - /j)Tjo(A)Tj '(/i) = T j0 >(^)T j0 (X)RoQ>(X - y) 



(5.10) 



where R(X) is the quantum R matrix associated to r 




(5.11) 
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approach to the question of integrable defects in classical field theories and, by no means, 
represents a complete picture of the story. 

Rather, it is a first step for future developments among which further study of the 
classical r matrix approach and quantization of the method are important. Let us mention 
also the contruction of other integrable defects allowing if possible reflection as well. If 
applicable, the quantization should then be related to existing quantum algebraic frame- 
works like the Reflection- Transmission algebras [6]. Finally, the complete setup of the 
direct and inverse part of the method for the actual construction of the solutions, espe- 
cially of soliton type, should shed new light on the results already obtained by the more 
direct approach of [HI [12] . 
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